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Thermal  energy  from  different  sources  can  be  converted  into  electricity  by  using  thermoelectric  modules. 
Controlling  the  temperature  uniformity  over  the  heat  spreader  of  the  thermoelectric  power  generation 
system  is  a  critical  issue  for  its  performance.  In  this  study,  we  compared  the  outputs  of  the  thermoelectric 
power  generation  system  which  has  uniform  temperature  distribution  with  that  without  uniform  tem¬ 
perature  distribution  over  the  heat  spreader.  The  model  in  which  thermoelectric  materials  are  tempera¬ 
ture  related  was  developed.  Finite  element  based  software  with  thermal-electrical  multi-physics  coupled 
was  employed  to  simulate  the  system.  In  addition  to  the  simulation,  we  performed  experimental  activ¬ 
ities.  In  the  simulation  and  experimental  activities,  electric  current  output,  power  output,  open  circuit 
voltage  and  the  efficiency  of  the  system  in  the  cases  of  uniform  and  non-uniform  temperature  distribu¬ 
tion  over  the  heat  spreader  are  evaluated.  The  simulation  and  experimental  results  show  that  the  uniform 
temperature  distribution  over  the  heat  spreader  gives  better  outputs  than  the  non-uniform  temperature 
distribution. 
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1.  Introduction 

The  direct  conversion  of  heat  energy  into  electricity,  or  the 
reverse,  by  a  semi-conductor  thermoelectric  power  generation  de¬ 
vices  is  related  to  electron  transport  phenomena,  and  the  interre¬ 
lated  Seebeck,  Petlier  and  Thomson  effects  [1,2].  Thermoelectric 
modules  are  devices  that  generate  voltage  and  electrical  power 
by  means  of  temperature  difference  [3].  They  are  of  great  interest 
in  energy  applications  due  to  their  well-known  merits  such  as  the 
absence  of  moving  components,  a  reduction  of  maintenance,  and 
an  increase  of  system  life.  Its  modularity  allows  for  application  in 
a  wide  scale  range  without  significant  losses  in  performance;  the 
absence  of  working  fluid  avoids  environmental  dangerous  leakage 
[4-8].  They  are  also  direct  current  (DC)  sources  and  have  capability 
of  recovering  the  huge  amount  of  low  grade  waste  heat  [7].  They 
comprise  p-  and  n-type  semiconductor  legs  connected  thermally 
in  parallel  and  electrically  in  series,  and  sandwiched  between 
two  ceramic  hot  and  cold  plates  [9-12]. 

Many  researchers  concerned  about  the  physical  properties  of 
thermoelectric  material  and  the  manufacturing  technique  of  ther¬ 
moelectric  modules.  In  addition  to  the  research  of  thermoelectric 
material  properties,  reasonable  thermal  system  design  and  optimi¬ 
zation  of  thermoelectric  generator  are  equally  important  for 
improving  the  power  generation  performance.  Rodriguez  et  al. 
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[13  developed  a  computational  model  to  simulate  the  thermal 
and  electrical  behavior  of  thermoelectric  generators.  Their  model 
solves  the  nonlinear  system  equations  of  the  thermoelectric  and 
heat  transfer.  Ahmet  and  Bekir  [14]  carried  out  a  theoretical  anal¬ 
ysis  of  thermoelectric  power  generator  and  they  formulated  the 
influence  of  thermoelectric  leg  geometry  on  the  device  efficiency 
and  power  generation.  Riffat  and  Ma  [15]  performed  the  geometry 
optimization  of  the  thermoelectric  modules  used  as  generator. 
Chen  et  al.  16]  and  Yu  et  al.  [17]  explored  the  influence  of  various 
parameters  on  the  performances  including  the  effect  of  multiple 
layers  of  modules.  Xuan  et  al.  [18  employed  a  phenomenological 
model  to  study  the  effects  of  internal  and/or  external  interface  lay¬ 
ers  on  the  thermoelectric  devices  performance.  Crane  et  al.  [11] 
investigated  the  characteristics  of  single-stage,  and  multi-element 
thermoelectric  generators  with  the  irreversibility  of  finite-rate 
heat-transfer,  Joulean  heat  inside  the  thermoelectric  device,  and 
the  heat  leak  through  the  thermoelectric  couple  legs.  Nuwayhid 
et  al.  [19]  analyzed  the  effect  of  the  finite-rate  heat  transfer  be¬ 
tween  the  thermoelectric  device  and  its  external  heat  reservoir 
on  the  performance  of  single-element,  and  single-stage  thermo- 
electric-generator. 

Some  researchers  used  thermoelectric  generators  in  their  model 
as  an  irreversible  heat  engine  [20].  Gao  and  Rowe  [21],  Rowe  and 
Min  [22]  and  Elena  et  al.  [23]  investigated  the  applications  of  ther¬ 
moelectric  power  generation  using  an  industrial  waste  heat  and 
low  temperature  waste  heat  to  generate  electric  power.  Some 
researchers  explored  the  integrated  thermoelectric  generators  in 


0196-8904/$  -  see  front  matter  ©  2013  Elsevier  Ltd.  All  rights  reserved. 
http://dx.doi.Org/10.1016/j.enconman.2013.08.004 


534 


B.T.  Admasu  et  al.) Energy  Conversion  and  Management  76  (2013)  533-540 


Nomenclature 

C 

specific  heat  capacity  (J  kg-1  K-1) 

Apn 

Ap  +  An 

/ 

electrical  scalar  potential  (V) 

P 

density  (kg  m-3) 

I 

electrical  current  (A) 

G 

electrical  conductivity  (S  m-1) 

j 

electrical  current  density  (A  m-2) 

T 

Thomson  coefficient  (V  K_1) 

N 

number  of  couple  legs 

P 

power  (W) 

Subscripts 

Q 

heat  generation  (W  m-3) 

av 

average 

R 

electrical  resistivity  (Q  m) 

c 

cold  junction 

Rpn 

Rp  +  Rn 

ce 

ceramic 

T 

temperature  (K) 

cu 

copper 

V 

open  circuit  voltage  (V) 

h 

hot  junction 

n 

n-type  module  leg 

Greek  symbols 

P 

p-type  module  leg 

a 

Seebeck  coefficient  (V  K  a 

Mpn 

o(p  an 

Abbreviations 

A 

change 

DC 

direct  current 

8 

dielectric  permittivity  (F  m-1) 

te 

thermoelectric 

n 

system  efficiency  (%) 

A 

thermal  conductivity  (W  m_1  K-1) 

various  micro-reactors  and  micro-combustors  [24,25].  Siddig  et  al. 
[26],  Chen  [27]  and  Scherrer  et  al.  [28]  performed  different  re¬ 
searches  on  solar  thermoelectric  power  generation  using  different 
system  configuration  and  arrangements. 

The  conversion  efficiency  of  thermoelectric  power  generation  is 
mainly  restricted  by  thermoelectric  materials.  Over  the  past  sev¬ 
eral  years,  a  number  of  high  power  performance  thermoelectric 
materials  have  been  developed  29,30],  and  some  of  them  are 
available  commercially  [31,32].  Sevan  et  al.  [33]  investigated  the 
characteristics  of  a  thermoelectric  generator  at  low  temperatures. 
Experimentally,  they  explored  the  temperature  dependency  of  See- 
beck  coefficient  and  electrical  conductivity  of  a  bismuth  telluride 
based  thermoelectric  generator. 

So  far,  there  are  so  many  researches  on  thermoelectric  power 
generating  system,  but  most  of  them  did  not  consider  the  effects 
of  temperature  non-uniformity  on  the  outputs  of  the  thermoelec¬ 
tric  power  generation  system.  The  purpose  of  the  present  study 
is  to  realize  this.  The  study  compares  the  outputs  of  thermoelectric 
system  when  the  input  temperature  is  uniformly  distributed  over 
the  heat  spreader  with  non-uniform  temperature  distribution  case. 
The  analysis  model  considers  the  temperature  dependency  of  all 
the  important  thermoelectric  material  properties.  The  second  and 
third  sections  of  this  study  introduce  the  governing  equations 
and  finite  element  modeling  of  thermoelectricity,  respectively. 
The  fourth  section  introduces  the  boundary  conditions.  The  fifth 
section  introduces  the  material  properties.  The  sixth  section  intro¬ 
duces  the  mesh  independency  of  the  thermoelectric  power  gener¬ 
ation  system.  The  seventh  section  introduces  the  results  and  the 
discussion  part.  The  final  part  presents  the  conclusion  part  of  this 
paper. 

2.  Governing  equations  of  thermoelectricity 

The  general  heat  conduction  and  continuity  equations  of  elec¬ 
tric  charge  for  the  thermoelectric  analysis  can  be  expressed  as 

[34-36]: 

pC^  +  V-m-J)-V-(W-VT)  =  q,  (1) 

V'  ([el'V|)  +  V'([<7]'[a]'V;r)  +  V'([<7]'V/)  =  0  (2) 


where  p ,  C,  T,J,f  and  q  are  the  density,  specific  heat  capacity,  abso¬ 
lute  temperature,  electric  current  density  vector,  electric  scalar 
potential,  and  heat  generation  rate  per  unit  volume,  respectively. 
[2]  is  the  thermal  conductivity  matrix,  [g]  is  the  electric  conductivity 
matrix,  [a]  is  the  Seebeck  coefficient  matrix,  and  [s]  is  the  dielectric 
permittivity  matrix. 

In  a  steady  state  thermoelectric  analysis,  thermal  and  electrical 
loads  do  not  vary  over  time.  In  the  present  steady-state  model, 
material  properties  of  all  components  are  considered  to  be  an  iso¬ 
tropic.  The  coupled  equations  of  thermoelectricity  can  be  written 
as  follows: 

V  •  (aocVT)  +  V  •  (aVf)  =  0  (3) 


V  •  (Taj)  -  V  •  (AVT)  =  q  (4) 

The  heat  generation  rate,  q  is  expressed  by  [37] 

q=J-  a VT  -  V  •  (AVT)  -hzJVT  (5) 

where  z  is  the  Thomson  coefficient  and  it  can  be  expressed  as 
follows 


a 

i  =  t3t 


(6) 


Keeping  the  non-ohmic  current  voltage  relation  in  both  legs  of 
the  devices  in  mind,  the  governing  equations  of  the  electrical  field 
inside  a  control  volume  under  steady  state  can  be  written  as  [4] 


(VV)p  =  -ccp(Tp)VTp-Rp(Tp)J 


(7) 


(W )n  = -ocn(Tn)VTn  -  Rn(Tn)J  (8) 

The  first  terms  in  the  right  side  of  Eqs.  (7)  and  (8)  are  the  Seebeck 
electromotive  force  (EMF)  increase  due  to  the  temperature  gradient, 
and  the  second  terms  are  the  voltage  drop  due  to  the  current  flow¬ 
ing  through  the  control  volume. 


3.  Finite  element  modeling 

Fig.  1  shows  the  finite  element  model  of  the  thermoelectric 
power  generation  system.  Here,  material  selections,  meshing  and 
boundary  conditions  for  the  model  were  set.  The  thermo-electric 
analysis  option  within  the  ANSYS  Workbench  was  used  to 
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Fig.  1.  Finite  element  model  of  thermoelectric  power  generation  system. 


represent  the  bismuth  telluride  as  thermoelectric  couple  legs 
material.  The  heat  spreader,  dielectric  ceramic  cover  and  the  heat 
sink,  which  are  the  other  components  of  the  thermoelectric  power 
generation  system,  included  in  the  modeled  as  thermal  masses 
with  isotropic  thermal  conductivity  and  electrical  resistivity. 

4.  Boundary  conditions 

The  study  set  different  boundary  conditions  to  simulate  the  re¬ 
search  model.  The  input  temperature  over  the  heat  spreader  was 
variable  while  the  heat  sink  temperature  was  an  average  ambient 
temperature  of  298  K.  The  free  surfaces  of  the  heat  spreader  and 
the  heat  sink  were  exposed  to  the  surrounding  air  and  heat  might 
be  removed  from  the  surfaces  due  to  natural  convection.  The  film 
coefficient  for  these  surfaces  assumed  to  be  25  W/(m2  K)  at  an 
average  ambient  temperature.  The  free  surfaces  of  the  thermoelec¬ 
tric  legs,  the  connectors  and  the  dielectric  ceramic  covers  were 
insulated  and  had  an  almost  zero  film  coefficient  that  is 
1  x  10-8  W/(m2  K)  at  an  average  ambient  temperature.  Initially,  a 
DC  voltage  of  0  V  (electrically  grounded)  and  0.0006  V  set  for  the 
low  potential  and  the  high  potential  boundary  conditions, 
respectively. 

5.  Material  properties 

The  selection  of  thermoelectric  materials  directly  affects  the 
performance  of  the  thermoelectric  device.  The  study  used  bismuth 
telluride  as  material  of  thermoelectric  module  legs.  The  properties 
of  this  material  such  as  the  Seebeck  coefficient,  thermal  conductiv¬ 
ity  and  electrical  conductivity  all  depend  on  temperature.  The 
physical  properties  of  the  thermoelectric  module  legs  provided 
by  Melcor  [38]  were  used  for  this  simulation  and  shown  as  follows 

0cpn  =  (22224  +  930.6Tav  -  0.9905T^V)  x  1(T9  (V/I<)  (9) 

Rpn  =  (5112  +  163.4T„V  +  0.6279T^V)  x  10“10  (Qm)  (10) 

V  =  (62,605  -  211 11  m  +  0.4131T^V)  x  10“4  (W/(mI<))  (11) 

where  TqV  =  ( Th  +  Tc)l 2  is  the  mean  temperature  of  the  hot  and  cold 
junctions. 

The  material  type  of  the  connectors,  heat  spreader  and  the  heat 
sink  for  this  simulation  was  copper.  The  physical  properties  of  cop¬ 
per  and  ceramic  at  an  average  temperature  of  400  K  are  shown  in 

Table  1. 

6.  Mesh  independency 

The  study  performed  a  mesh  sensitivity  analysis  for  the  thermo¬ 
electric  power  generation  system  to  check  its  mesh  independency. 


Table  1 

Physical  properties  of  copper  and  ceramic. 


Physical  properties 

Units 

Copper 

Ceramic 

References 

X 

W/(m  K) 

398 

25.12 

[39] 

R 

Q  m 

i.67  x  nr8 

5.03  x  1010 

[39] 

The  default  method  for  meshing  that  is  rectangular  mesh,  which  is 
shown  in  Fig.  2.  The  study  prepared  three  models  to  check  the 
mesh  sensitivity  and  the  mesh  size  was  progressively  refined. 
The  first  model  had  37,760  nodes  and  4959  elements  with  element 
size  of  9  x  10-4  m,  and  yields  a  current  output  of  5.6077  A,  the  sec¬ 
ond  model  had  41,136  nodes  and  5427  elements  with  element  size 
of  8.5  x  10-4  m,  and  yields  a  current  output  of  5.6034  A,  and  the 
third  model  had  44,542  nodes  and  5890  elements  with  a  mesh  size 
of  8  x  10-4  m,  and  yields  a  current  output  of  5.6017  A.  Finally,  the 
study  confirmed  that  as  the  mesh  size  refines  more,  the  current 
output  does  not  have  obvious  change. 

7.  Results  and  discussion 

7 A.  Numerical  simulation 

The  performance  of  thermoelectric  power  generation  system 
might  be  negatively  affected  by  different  factors,  such  as  losses 
due  to  thermal  resistance  of  the  components,  losses  due  to  thermal 
contact  and  electrical  contact  resistances  at  the  interfaces  between 
the  components,  and  temperature  non-uniformity  on  the  heat 
spreader.  All  the  above  losses  are  due  to  the  geometry  of  the  ther¬ 
moelectric  power  generation  system.  In  the  present  paper,  we  fo¬ 
cused  on  the  effect  of  temperature  non-uniformity.  As  seen  from 
the  finite  element  model  in  Fig.  1,  the  heat  spreader  was  divided 
into  four  partitions.  In  the  first  case  that  is  the  uniform  tempera¬ 
ture  distribution,  the  study  used  the  same  temperature  values  for 
the  four  partitions  and  the  outputs  were  evaluated  for  different 
temperature  ranges.  In  the  second  case  that  is  the  non-uniform 
temperature  distribution,  the  study  used  different  temperature 
values  at  the  four  partitions  but  an  equal-average  temperature 
with  that  of  the  case  of  the  uniform  temperature  distribution  at 
each  step.  The  other  inputs  were  the  geometry  data,  the  boundary 
conditions  and  temperature  dependent  thermoelectric  material 
properties  such  as  the  Seebeck  coefficient,  thermal  conductivity 
and  electrical  resistivity. 

Figs.  3  and  4  show  the  simulation  results  of  the  temperature 
contour  in  the  cases  of  uniform  and  non-uniform  temperature  dis¬ 
tribution  over  the  heat  spreader.  The  sample  temperature  data  ta¬ 
ken  to  simulate  in  the  case  of  uniform  temperature  distribution 
were  360.75  K,  370.75  K,  380.75  K,  390.75  K,  400.75  K,  410.75  K, 
420.75  K  and  430.75  K  and  the  temperature  data  in  the  case  of 
non-uniform  temperature  distribution  were  (435  K,  360  K,  350  K 
and  298  K),  (435  K,  380  K,  370  K  and  298  K),  (435  K,  400  K,  390  K 

and  298  K),  (435  K,  420  K,  410  K  and  298  K),  (435  K,  440  K,  430  K 

and  298  K),  (435  K,  460  K,  450  K  and  298  K),  (435  K,  480  K,  470  K 

and  298  K)  and  (435  K,  500  K,  490  K  and  298  K).  The  average  tem¬ 

perature  of  the  non-uniform  temperature  distribution  at  each  step 
was  equal  with  that  of  the  uniform  temperature  distribution.  For 
example,  the  first  input  temperature  data  in  the  case  of  uniform 
temperature  distribution  that  is  360.75  K  is  equal  with  the  average 
temperature  of  the  first  input  temperature  data  of  the  non-uniform 
temperature  distribution  that  is  (435  K,  360  K,  350  K  and  298  K). 

The  hot  junctions  of  all  the  thermoelectric  module  couple  legs 
are  inter-connected  by  thermal  conductive  copper  strips.  When 
we  see  the  assembly  of  the  system,  all  the  conductive  copper  strips 
are  inter-connected  with  the  upper  and  lower  dielectric  ceramic 
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Fig.  2.  Mesh  model  of  thermoelectric  power  generation  system  with  rectangular  elements. 


370.78  362.69  354.6  346.51  338.42  330.33  322.25  314.16  306.07  297.98 


Fig.  3.  Simulation  result  of  temperature  contour  in  the  case  of  uniform  input  temperature  of  370.75  K. 


Fig.  4.  Simulation  result  of  temperature  contour  in  the  case  of  non-uniform  input  temperature  of  (435  K,  380  K,  370  K  and  298  K). 
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Fig.  5.  Voltage  output  versus  the  temperature  change  between  the  junctions  in  the  Fig.  8.  Power  output  and  system  efficiency  versus  the  temperature  change  between 

cases  of  uniform  and  non-uniform  temperature  distribution  over  the  heat  spreader.  the  junctions  in  the  cases  of  uniform  and  non-uniform  temperature  distribution 

over  the  heat  spreader. 


Fig.  6.  Electric  current  output  versus  the  temperature  change  between  the 
junctions  in  the  cases  of  uniform  and  non-uniform  temperature  distribution  over 
the  heat  spreader. 


Fig.  7.  Electric  current  density  versus  the  temperature  change  between  the 
junctions  in  the  cases  of  uniform  and  non-uniform  temperature  distribution  over 
the  heat  spreader. 


cover  which  have  direct  contact  with  the  heat  spreader  and  heat 
sink,  respectively.  One  can  see  from  the  temperature  contour  of 
Fig.  4a  and  b,  in  the  case  of  non-uniform  temperature  distribution, 
on  a  single  system,  the  hot  junction  of  all  the  module  legs  have  dif¬ 
ferent  temperature  values.  As  seen  in  Fig.  4b,  even  on  a  single  mod¬ 
ule  leg,  at  different  coordinates,  the  temperature  distribution  is 
non-uniform.  As  an  example,  on  a  single  module  leg,  the  hot  junc¬ 
tion  temperature  values  at  different  coordinates  are  431.06  K  and 
425.84  K.  One  can  see  from  Fig.  4  that  some  parts  of  the  heat  sprea¬ 
der  have  lower  temperature  values  than  some  parts  of  the  module 
legs  hot  junction  temperature  values.  So,  due  to  the  second  law  of 
thermodynamics,  heat  might  be  flowed  from  the  hot  junction  of 
the  module  legs  to  the  heat  spreader  which  is  in  the  reverse  direc¬ 
tion.  If  the  temperature  distribution  on  the  heat  spreader  is  uni¬ 
form,  as  seen  in  Fig.  3a  and  b,  heat  will  be  flowed  only  from  the 
heat  spreader  to  the  hot  junction  of  the  module  legs. 

The  junctions’  temperature  cannot  be  determined  simply  from 
the  input  temperature  on  the  heat  spreader  due  to  the  thermoelec¬ 
tric  effect  inside  the  thermoelectric  module  legs.  Those  thermo¬ 
electric  effects  are,  the  internal  irreversibility  caused  by  Joulean 
electrical  resistivity  loss,  heat  conduction  loss  through  the  semi¬ 
conductor  legs  between  the  hot  and  cold  junctions,  Peltier  effect 
at  the  hot  and  cold  junctions,  and  a  reversible  absorption  or  liber¬ 
ation  of  heat  in  a  homogenous  temperature  and  electrical  gradient 
for  temperature  dependent  thermoelectric  materials  due  to  Thom¬ 
son  effect.  So,  reading  the  temperature  from  the  simulation  result 
at  the  junctions  was  indispensable. 

To  simulate  the  system  for  an  open  circuit  voltage,  low  potential 
boundary  condition  was  set  on  the  external  face  of  conductive  cop¬ 
per  strip  that  connects  the  external  electrical  load  and  the  n-leg, 
and  the  high  potential  boundary  condition  was  set  on  the  external 
face  of  the  conductive  copper  strip  that  connects  the  external  elec¬ 
trical  load  and  the  p-leg. 

Fig.  5  shows  the  thermoelectric  power  generation  system  volt¬ 
age  output  versus  the  change  in  temperature  between  the  junc¬ 
tions  in  the  cases  of  uniform  and  non-uniform  temperature 
distribution  over  the  heat  spreader.  It  can  be  seen  that  the  voltage 
output  is  linear  and  increasing  function  of  the  change  injunctions’ 
temperature.  The  linear  relationship  is  approximate  because  there 
are  factors  like  temperature  dependent  thermoelectric  material 
properties  and  radiative  heat  transfer  which  caused  non-linearity. 
As  seen  from  Fig.  5,  when  the  temperature  change  between  the  hot 
and  cold  junctions  increases,  the  voltage  outputs  in  both  cases 
increase  but  the  increment  in  the  case  of  uniform  temperature 


538 


B.T.  Admasu  et  al.) Energy  Conversion  and  Management  76  (2013)  533-540 


(c)  (d) 


Fig.  9.  Experimental  components  and  setup  (a)  the  heater  covered  all  top  surface  of  the  module,  (b)  the  heater  covered  some  parts  of  the  module,  (c)  and  (d)  experiment  set 
up  in  case  (a)  and  case  (b),  respectively. 


distribution  is  relatively  higher.  We  can  see  from  Fig.  5  that  when 
the  temperature  change  is  90  K,  the  voltage  outputs  in  the  cases  of 
uniform  and  non-uniform  temperature  distribution  over  the  heat 
spreader  are  0.567  V  and  0.55  V,  respectively  and  the  voltage 
difference  between  the  two  cases  is  0.017  V.  When  the  tempera¬ 
ture  change  is  128  K,  the  voltage  output  in  the  cases  of  uniform 
and  non-uniform  temperature  distribution  are  0.84  V  and  0.81  V, 
respectively  and  the  voltage  difference  between  the  two  cases  is 
0.03  V. 

The  electric  current  generated  has  direct  relationship  with  the 
open  circuit  voltage.  Based  on  the  Seebeck  effect,  and  the  temper¬ 
ature  difference  between  the  hot  and  cold  junctions,  causes  an 
electric  current  to  flow  in  the  circuit.  Fig.  6  shows  the  electric  cur¬ 
rent  output  versus  the  temperature  change  between  the  junctions 
in  the  cases  of  uniform  and  non-uniform  temperature  distribution 
over  the  heat  spreader.  As  seen  from  the  figure,  when  the  temper¬ 
ature  change  between  the  hot  and  cold  junctions  increases,  the 
electrons  on  the  hotter  side  vibrates  more  vigorously  so  they  tend 
to  move  towards  the  colder  side  where  the  electrons  are  moving 
more  slowly  and  the  current  generation  will  increase.  From 
Fig.  6,  one  can  see  that  the  electric  current  is  parabolic-like  and 
increasing  function  of  the  temperature  change  between  the  hot 
and  cold  junctions  in  both  uniform  and  non-uniform  temperature 
cases.  Actually,  the  current  output  and  temperature  have  linear 
relationships  but  there  are  factors  like  temperature  dependence 
of  thermoelectric  material  properties  and  radiative  heat  transfer 
to  cause  the  non-linearity.  As  mentioned  in  the  previous  section, 
due  to  the  second  law  of  thermodynamics,  the  current  output  is 
relatively  higher  in  the  case  of  uniform  temperature  distribution 
over  the  heat  spreader  than  an  equal  average  non-uniform  temper¬ 
ature  distribution.  Fig.  7  shows  the  electric  current  density  versus 
the  change  in  temperature  between  the  junctions  in  the  cases  of 
uniform  and  non-uniform  temperature  distribution  over  the  heat 
spreader.  It  has  the  same  profile  with  that  of  the  electric  current 
output  because  all  the  parameters  influence  the  current  density. 
As  the  temperature  gradient  increases,  the  electric  current  density 


increases  parabolic-like.  As  seen  from  the  figure,  if  we  keep  the 
temperature  distribution  uniform  over  the  heat  spreader,  we  can 
get  better  particle  flow. 

The  electric  power  output  can  be  evaluated  from  the  heat  input 
at  the  hot  junction  and  the  heat  dissipated  from  the  thermoelectric 
module  through  the  heat  sink  and  the  efficiency  of  the  system  can 
be  evaluated  from  the  ratio  of  electric  power  output  to  the  heat  in¬ 
put  at  the  hot  junction  of  the  thermoelectric  module.  Fig.  8  shows 
the  power  output  and  system  efficiency  versus  the  temperature 
change  between  the  junctions  in  the  cases  of  uniform  and  non-uni- 
form  temperature  distribution  over  the  heat  spreader.  It  can  be 
seen  that  the  power  output  and  efficiency  are  parabolic-like  and 
increasing  function  of  the  change  in  temperature  between  the 
hot  and  cold  junctions  in  both  cases  that  is  the  uniform  and  non- 
uniform  temperature  distribution  over  the  heat  spreader.  The 
power  output  and  system  efficiency  increment  in  the  case  of  uni¬ 
form  temperature  distribution  are  relatively  higher.  One  can  see 
from  Fig.  8  that  when  the  change  in  temperature  is  80.5  K,  the 
power  output  and  system  efficiency  are  1.562  W  and  3.81%  in 
the  case  of  uniform  temperature  distribution  and  1.425  W  and 
3.71%  in  the  case  of  non-uniform  temperature  distribution,  respec¬ 
tively.  The  differences  in  power  output  and  system  efficiency  at 
this  point  between  the  two  cases  are  0.137  W  and  0.1%,  respec¬ 
tively.  When  the  change  in  temperature  is  129  K,  the  power  output 
and  system  efficiency  are  4.68  W  and  5.77%  in  the  case  of  uniform 
temperature  distribution  and  4.032  W  and  5.62%  in  the  case  of 
non-uniform  temperature  distribution,  respectively.  The  differ¬ 
ences  in  power  output  and  system  efficiency  at  this  point  between 
the  two  cases  are  0.648  W  and  0.15%,  respectively. 

7.2.  Experimental  validation 

The  commercial  thermoelectric  module  TEC1 -12706  type  and 
the  resistance  wire  silicon  rubber  heater  TY-SRH  were  tightly 
clamped  between  the  upper  frame  and  the  circulating  water  reser¬ 
voir.  The  silicon  rubber  heater  glued  to  the  module  using  epoxy 
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the  data  acquisition  multimeter,  “K”  type  thermocouples  were 
used  to  measure  the  temperatures  of  heat  sink,  the  hot  and  cold 
surfaces  of  the  thermoelectric  module,  the  heater  and  the  ambient. 

The  Circulating  water  used  for  heat  sink  had  a  maximum  flow 
rate  of  8.3  L  min-1  and  a  maximum  head  of  0.8  m.  The  study  used 
two  different  sample  sizes  of  silicon  rubber  heaters.  The  first  sam¬ 
ple  covered  the  whole  upper  surface  of  the  thermoelectric  module 
as  seen  from  Fig.  9a  which  has  the  same  size  with  the  thermoelec¬ 
tric  module,  40  mm  by  40  mm  and  delivers  uniform  temperature 
distribution.  As  seen  from  Fig.  9b,  the  second  sample  size  was 
15  mm  by  40  mm  and  covered  and  heat  up  some  part  of  the  upper 
surface  of  the  thermoelectric  module.  When  the  experiment  per¬ 
formed,  we  used  the  insulating  material  (i.e.  asbestos)  to  cover 
the  surroundings  of  the  heater  and  the  uncovered  parts  of  the 
module. 

In  the  experiment,  the  parameters  for  the  two  cases  were  the 
same  only  the  coverage  of  the  heater.  When  we  performed  the 
experiment,  we  connected  an  external  electrical  load  of  1.5  Q. 
Fig.  10  shows  the  experimental  results  on  the  effect  of  temperature 
non-uniformity  over  the  heat  spreader.  As  seen  from  the  figure,  if 
we  keep  the  temperature  distribution  uniform  over  the  hot  surface 
of  the  thermoelectric  module,  we  can  get  better  output. 

8.  Conclusion 

In  this  paper,  we  formulated  the  governing  equations  of  ther¬ 
moelectricity  and  conducted  finite  element  modeling  and  analysis 
to  compare  the  effect  of  uniform  and  an  equal  average  non-uni- 
form  temperature  distribution  on  the  output  of  thermoelectric 
power  generation  system.  The  experimental  and  simulation  results 
show  that  maintaining  the  temperature  distribution  uniform  all 
over  the  top  surface  of  the  heat  spreader  of  the  thermoelectric 
power  generation  system  delivers  better  outputs. 
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